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The thermodynamics of the spin-5* anisotropic quantum XXZ chain with arbitrary value of S 
and unitary norm, in the high-temperature regime, is reported. The single-ion anisotropy term and 
,— _l ■ the interaction with an external magnetic field in the z-direction are taken into account. We obtain, 

D ' for arbitrary value of 5*, the /3-expansion of the Helmholtz free energy of the model up to order 

5_j , /3 6 and show that it actually depends on S ^+V) • ^ s classical limit is obtained by simply taking 

■ S — > oo. At h = and D — 0, our high temperature expansion of the classical model coincides 
with Joyce's exact solution [TUl . We study, in the high temperature region, some thermodynamic 
quantities such as the specific heat and the magnetic susceptibility as functions of spin and verify 
for which values of S those thermodynamic functions behave classically. Their finite temperature 
behavior is inferred from interpolation of their high- and low-temperature behavior, and shown to 

"-^J , be in good agreement with numerical results. The finite temperature behavior is shown for higher 

C"j ■ values of spin. 

o 

Q ; PACS numbers: 05, 05.30.Ch, 75.10.Jm, 75.10.Hk, 65. 

I. INTRODUCTION 

i> ; 

One-dimensional and quasi-one-dimensional quantum chain models have been intensively investigated for a twofold 
reason. On one hand, for the algebraic properties of spin-1/2 chain models {XXZ and generalized ladder models) 
which, being integrable models, can be solved exactly by the powerful Bethe ansatz technique (see Refs. pj and 0], 
ly— j ■ respectively). On the other hand, for the existence of quasi one-dimensional magnetic systems that exhibit nearly 
ideal one-dimensional behavior for some interval of temperature. The thermodynamic properties of those magnetic 
systems have been described by spin-5 1 XXZ models for various values of spin, mixed-spin models, ladder models, 
etc. In particular, there are materials well described by higher spin values of the XXZ models such as CSVCI3 and 
CsVBrg {S = 3/2)|l|ll3, (Ci H 8 N 2 )MnCl 3 (S = 2)0 and (CD 3 ) 4 NMnCl 3 (S = 5/2)0 @. Motivated by the 
discovery of these materials, the magnetic and thermodynamic properties of the ferro- and antiferromagnetic isotropic 
spin-5 XXZ model for higher values of S have been studied numerically and analytically (e.g., the high temperature 
expansions; see Ref. and references therein). 

The classical limit (S — > 00) of the XXZ model is exactly soluble and was originally solved by Fisher^cji m t ne 
isotropic regime, for a null external magnetic field and no single-ion anisotropy term. Its anysotropic case has been 
solved by Joyce [Til ], also in the absence of an external magnetic field and without a single- ion anisotropy term, by 
writing the Helmholtz free energy of the classical XXZ model as an integral equation whose solution is the spheroidal 
wave function. Unfortunately, in the presence of an external magnetic field and/or the single- ion anisotropy parameter, 
the solution of the integral equation cannot be reduced to any known function. 

Recently, Fukushima et al.]l^ obtained the specific heat and the magnetic susceptibility per site of a ferromagnetic 
mixed-spin model, with two kinds of spins, s and S, arranged alternatively and coupled by a Heisenberg-type nearest- 
neighbor exchange for arbitrary values of s and S at the isotropic point (A = 1 in cq.Q) and in the absence of external 
magnetic fields. They obtained numerical results through the exact diagonalization method and the analytical high 
temperature expansions. For s = S their model becomes a one-dimensional spin-S 1 XXZ model for arbitrary value of 
S. Their high temperature expansion of the specific heat per site, in the absence of an external magnetic field, goes 
up to order (/3J) 11 and their expansion of the magnetic susceptibility, also for a vanishing magnetic field, goes up to 
(PJ) 7 with a single-ion anisotropy term in one of the spins. In Ref. [jj we extended the results of Ref. [T3I, for spin 
values up to S = 4, including a free parameter of anisotropy in the z-direction, a single-ion anisotropy term and a 
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non-zero external magnetic field in the /3-expansions of the Helmholtz free energies up to order (J(3) 6 . 

Following Ref. 0, we also consider the Hamiltonian of the anisotropic spin-5 XXZ with a single- ion anisotropy 
term and in the presence of an external magnetic field, 



JV 

H = £ J' (Si, S i+1 ) A - h'St + D'iSff. (1) 

i=l 

We use the notation: (S,-, S i+1 ) A = SfSf +1 + SfSf +1 + A5?5? +1 . Here, Sf, Sf and Sf stand for the spin-5 matrices 
in the i-th site of the chain and norm */ 5(5 + 1); N is the number of sites in the periodic chain; J' is the exchange 
integral; A is the anisotropy constant in the z-direction; h! is the external magnetic field in the z-axis and D' is the 
single-ion anisotropy parameter. 

Even in one spatial dimension, the solution of the quantum chain model Q is more complex if we consider a finite 
arbitrary spin-5 (with 5 ^ 1/2) once the XX Z model becomes a non-integrable model and cannot be solved by the 
thermodynamic Bethe ansatz technique. 

Numerical approaches for higher spin values becomes more involved since there are more degrees of freedom to be 
handled. Certainly, the high temperature expansions are easily manipulated by symbolic computer languages and can 
be used as a reliable check for numerical calculations in the high temperature regime. 

The aim of the present communication is to extend the results of Ref. .9] to arbitrary values of the spin 5 in 
hamiltonian (JIJ . In doing so, we are including in the results of Ref. [l2j the effects of anisotropies in the z-direction 
of the XX Z model in the presence of an external magnetic field. Having a high temperature expansion for the model, 
for any 5, we can derive its classical limit in this regime of temperature. We want to re-obtain the classical results of 
the Heisenberg model by considering the quantum nature of the spin variable, in opposition to the known results in 
the literature where this limit is derived from a chain of classical spins. From the high temperature expansion of the 
Helmholtz free energy, we can also verify for which values of 5 the anisotropic quantum XX Z modeL in the presence 
of an external magnetic field, is well described by its classical version, extending the results of Ref. Q • In the present 
article, comparison of thermodynamic functions for different values of 5 is performed at the same (high) T, whereas 
in Ref. [!j , they were plotted as functions of the "rescaled" temperature T = ■ 

This paper is organized as follows: in §11 we use the method developed in Ref. |l3| and an interpolation technique 
to obtain the high-temperature series expansion of the Helmholtz free energy for the quantum spin-s XXZ chain (s 
with unitary norm) up to order f3 6 . In order to check our analytic /3-expansion, in the first subsection of §11 our results 
for the classical limit (5 — ► oo) are compared to the well-known isotropic^^ an d anisotropic 0| classical Heisenberg 
chains. In the second subsection of §11 our results for finite values of 5 are explored to include the effect of the 
anisotropy parameter A, the presence of D' and b! in the Hamiltonian Q and the dependence of thermodynamic 
functions on 5. In §111 we use Pade representatives to enhance our high-temperature expansions for the spin-s XXZ 
model, extending them to lower temperatures. We do so by taking into account the known behavior at T ~ 0. For the 
specific heat we use the method presented in Ref. , whereas for the magnetic susceptibility we apply the Dlog-Pade 
approximant (antiferromagnetic case) and the two-point Pade approximant ^(J (ferromagnetic case). In §IV we 
present our conclusions. In appendix A we show that the trace of powers of S z can be written as an expansion in the 
square of the norm of the spin at each site. Finally, in appendix B we present the high-temperature expansion, up to 
order /3 6 , of the Helmholtz free energy of the spin-s XXZ model (|s| = 1), with a single- ion anisotropy term and in 
the presence of an external magnetic field. 



II. THE HIGH TEMPERATURE BEHAVIOR OF THE QUANTUM SPIN-s XXZ CHAIN 

We intend to study the thermodynamics of the Hamiltonian for arbitrary spin, including its classical limit (5 — > 
oo). In this limit, the thermodynamic functions diverge and to render them finite we define a rescaled spin operator 
s = S/y S(S + 1). This rescaled spin operator has unitary norm for all values of S. Rewriting the Hamiltonian Q 
in terms of s and redefining the parameters J = S(S + 1) J', h = y/S(S + l)h' and D = S(S + l)D', we obtain 

N 

H = ^J( Si , Sm ) A -K + D(sl) 2 . (2) 

i=l 

This also describes the dynamics of the XXZ model of a spin of unitary norm and (25 + 1) z-components. We 
continue to use the notation: (sj,Sj+i)A = s f s f+i + s i s ^+i + Asfsf +1 . Here, sf, sf and sf stand for the spin-s 
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(s = 1/2, 1, 3/2, • ■ ■ ) rescaled matrices in the i-th site of the chain; N is the number of sites in the periodic chain; and 
A is the anisotropy constant in the z-direction. 

In order to evaluate the Helmholtz free energy for arbitrary values of spin, we first calculate its expansion for a 
set of spin values using the Hamiltonian (J2J. We apply the method of Ref.pjl to calculate those /3-expansions for 
arbitrary values of J, A, h and D. We presented elsewhere|9j the Helmholtz free energy of Hamiltonian for a 
number of spin values (semi-integer and integer) up to S — 4. We use the interpolation method to extend our results 
for fixed values of S to arbitrary spin value. For the sake of obtaining the /3-expansion of any spin-5 quantum chain 
up to order /3 6 , having the expansions presented in Ref. 9] is not enough. 

In order to derive the expansion of the Helmholtz free energy of the hamiltonian (JIJ up to order f3 n (for a survey of 
the method (l3|. cf. section 2 of Ref. 0), we have to calculate normalized traces, at each site, of the products of m, 
matrices from the set {S x , S y , S z }, so that < 2(k + 1) and k ranges over all possible powers of /3, k = 1, 2, • • • , n. 
Due to the commutation relation [S x , S v ] = iS z and to the fact that the products S X S V and S V S X can be written in 
terms of S z , (S z ) 2 and the identity matrix 1, then the normalized trace at each site is reduced to the calculation of 
traces of (S z ) 1 , with I e N. 

In ea. ljAlj! we show that the tr((S z ) 1 ), for even values of I, is a polynomial of degree (/ + 1) in S. Therefore the 
polynomial in S of highest degree at order /3 n is of order 2(n + 1). In order to carry out the interpolation of the 
expansion up to order /3 6 , we calculated the /3-expansion of the Helmholtz free energy from the spin S — 1/2 up to 
5 = 7. Later, by inspection of the expansion in S and in /3, we verified that the expansion in S can be rewritten in 
terms of S{S + 1), that is, in terms of the eigenvalue of a constant of motion (the square of the norm of the spin-S at 
each site). In appendix 1X1 we show, for arbitrary value of spin, that tr[(Sf) ], I = 1, 2 • • • , is a polynomial of degree 
(/ — 1) in the parameter S(S + 1). 

The thermodynamic properties of the quantum spin-5 1 XX Z chain, in the high temperature region, can be obtained 
from its Helmholtz free energy. We present in appendix B the coefficients of the Helmholtz free energy W s (/3) derived 
from the Hamiltonian (J2J), up to order /3 6 , for any value of S(S = 1/2, 1, 3/2, • • • ). In particular, the result of appendix 
B is also valid at the isotropic point of the Hamiltonian (J2J) with D = and h — 0, when the model has rotational 
symmetry. The /3-expansion of the Helmholtz free energy of Hamiltonian (JIJ is easily obtained from result (|B3|I (see 
appendix [B|. Fukushima et aZ.|l2| obtained the /3-expansion of the specific heat and magnetic susceptibility per site 
of the latter model for arbitrary spin at the isotropic point and in the absence of an external magnetic field. Our 
results agree with theirs, up to order /3 6 . 

The range of validity in (3 of our high-temperature expansions obviously depends on the particular values given to 
the parameters of @. The thermodynamic properties of the chain depend on the sign of the product J A. We let 
J > and let A refer to either the ferromagnetic (A < 0) or the antiferromagnetic (A > 0) phases. 



A. Comparison with known classical limits (S — > oo) 

An interesting check to our expansion (|B3|I is to verify if it recovers known classical limits (S — > oo or y — > 0, where 
y = S (g +1 ) ) of the XX Z chain. The classical limit of the XX Z chain of the Hamiltonian J2J with h — and D = 

has been solved by Joyce [ijj; this model (including the anisotropic constant A in the Hamiltonian J2J) is exactly 
solvable and its Helmholtz free energy is given by the radial component of the spheroidal wave function 



Woo(m;/?) = ~0 in (-Ro,o(-«/3^/sinh fi, cosh /x)), (3) 

where /i = tanh -1 ^) and Rq,o is the first radial spheroidal wave function|17l]. For |A| = 1, eq. 10 becomes the 
well-known result obtained by Fisher^], whose explicit expression is Woo(±§ ,0) — — 4 ln[sinh(/?J)/(/3J)]. 

The classical limit of the /3-expansion (|B3|I is obtained by letting y = 0. In order to have a shorter /3-expansion of 
the Helmholtz free energy at y = 0, we take D = 0, and obtain 
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where h = h/J. 

From Q we see that in the absence of the external magnetic field (h — 0) the Helmholtz free energy is an even 
function of the anisotropic parameter A. Under this condition, the classical limit of the XXZ model does not 
distinguish between the ferromagnetic and the antiferromagnetic phases in the region of high temperatures. 




FIG. 1: The Helmholtz free energy Woo (for S — » oo) per unit of kT/J as a function of kT/J for the classical XXZ model in 
the absence of the single-ion anisotropy term (D = 0) and an external magnetic field (h = 0). In (a) we plot the anisotropic 
cases with A = ±0.5 and in (b) the isotropic cases (A = ±1). The solid line refers to the numerical solution of 10 and the 
dashed line to the expansion 



Fig. (JJJi) shows a comparison of our /3-expansion of the Helmholtz free energy to the numerical solution of eq.JSJ), 
for h = 0, D = and A = ±0.5. They are in good agreement in the high temperature region; the percental error 
between them at kT/J — 0.2 is less than 2.8%. Fig. QJ>) performs the same comparison, but at the isotropic point 
(A = ±1). We also have good agreement, but for higher temperatures; percental error is about 3% for kT/J = 0.4. 

In conclusion, the classical limit of the expansion (|jj3|) (for y = 0) coincides with Joyce's solution © of the 
Helmholtz free energy of the Hamiltonian @| in a larger interval of temperature than that which we call "high 
temperature region" . 
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B. The thermodynamics of the spin-s XX Z model 



The greater the (finite) spin, the greater the number of degrees of freedom to be handled, thus the more involved 
it gets to compute thermodynamic properties of the quantum spin-S model. In contrast, for infinite values of spin, 
the calculation of those properties for the classical version of the model is simpler, since we can treat the spin as a 
vector that rotates continuously. It is interesting to know for which values of S the quantum spin model can be well 
approximated by its classical version, at least in the high temperature region. 

In Ref. J|j we compared various thermodynamic functions of the chain model at the Heisenberg point (A = ±1), 
in the absence of the single-ion anisotropy term and external magnetic fields. We verified, in the high temperature 
region, that thermodynamic functions like the magnetic susceptibility per site have a classical behavior, for S > 2. 
Our intention in this section is to extend the results of Ref. |2( in order to include in this comparison (quantum x 
classical behavior) of spin models with anisotropics in the z-direction (A ^ ±1 and D ^ 0) and in the presence of an 
external magnetic field. This analysis can be realized from the expansion i|B3|l of W s ((3)- The referred thermodynamic 
functions are all derived from the Hamiltonian J5J. 

Let C s {(3) and Cs{[3) be the specific heat per site for the spin-s model (with unitary norm) and for the spin-S* 
chain (with norm equal to \J S(S + 1)), respectively. The specific heat per site is calculated from a derivative of the 

" (|3Wi 



Helmholtz free energy (C(/3) = -f3 2 d °3 (/3)) ). From the results JB2j| we obtain, 



(3 



C s (J,A,h,D;P) = C s U A, V ' S{S + l)h, D; M + - j (5a) 

= ^s^' A '7lWs^ ;/3 )- (5b) 

These results show that Cs is a homogeneous function of zero degree for all temperatures. The (3- expansion of the 
specific heat derived from eq. (|B3|) satisfies eqs.© identically. 

From the high temperature expansion of C s we verify that C s = —( j^sTs+T) ~ Ts — I ~~ TI lf)<^ 2 + ^(/3 3 )- in 

the high temperature region, the XX Z model presents a tail of the Schottky peak^^] (Csch oc /3 2 ), for all values of 
S. 

Fig. |U shows the specific heat per site as a function of J/3 for distinct values of S, including the classical limit 
(S — > oo) of the Hamiltonian J2J) and the relative percental error of this function for various values of S in relation to 
the classical specific heat per site. We take h/J = 0.5 and D/J = 0.7. Fig. J2t) pictures the antiferromagnetic case 
(A = 0.3) and (J2J;) the ferromagnetic case (A = —0.3). Figs. (J2J)) and |J2Ji) show the relative percental error, i.e., 

x 100%, of Figs. an d {30, respectively, where is the classical limit of the specific heat per site. In 



both cases, even the S = 2 model behaves classically (within an error smaller than 2%) up to J/3 = 0.3. Within this 
range of error, we can also neglect the quantum nature of the spin model with S > 5 up to J/3 ~ 1. 

Let Xs(P) and Xs(P) b e the magnetic susceptibilities per site of the spin-s model (with unitary norm) and of the 
spin-5 model (with norm equal to yj S(S + 1)), respectively. The relation between the magnetic susceptibilty per site 

>2 w( 

dh 2 



and the Helmholtz free energy is 

Xifi) = ~ 9 7v ] ■ From the results (JHSj we obtain the relation between \s and xs, 



Xs(J,A } h,D;[3) = X S U A, y/S(S + l)h, D; ^ j (6a) 

' ^fw^.A,^^,-^-;^. (6b) 



S(S + 1) M \S(S + 1)' ' ^5(5 + 1)' 5(5 + 1)' 

Equating the r.h.s. of eqs. I|6a|) and (|6b|l we obtain that xs is a homogenous function of degree 1 for all values of 
temperature. Its /3-expansion satisfies this property. 

Fig. 01 shows the magnetic susceptibility versus J/3 for various finite values of S and the classical chain model 
(S — > oo). Fig. JsB 3 ) shows the relative percentage error of the magnetic susceptibilty with respect to the classical 
curve for S = 1,3/2,4 and 6. In Fig. OJwe have A = 1, h/J = 0.3 and D/J = 0.5. From Fig. 03) we verify 
that the magnetic susceptibility of the spin-3/2 can be approximated by the classical result up to J (3 ~ 0.5 with an 
error smaller than 2%. The classical magnetic susceptibility curve is a good approximation of the spin-2 model up to 
J 13 ~ 0.7. 
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FIG. 2: Figs, (a) and (c) show the specific heat C a (J(3); (b) and (d) show the relative percental errors between C a (/3) and 
Coo(/3) (the classical limit) in (a) and (c), respectively. For all plots we have h/J = 0.5 and D/J — 0.7. In (a) we let A = 0.3 
and in (c) we have A = —0.3. In (a) and (c), line styles correspond to different values of spin, as follows: S = 1/2 (dashed and 
double-dotted lines), S — 1 (dotted line), S — 2 (dashed line), 5 = 3 (dashed and dotted lines) and S — > oo (solid line). In (6) 
and (d) we have S = 2 (solid line), S = 3 (dotted line), S = 5 (dashed line) and 5 = 10 (dashed and dotted line). 



Fig. Q]shows the magnetization of the XXZ model of the spin with unitary norm versus h/J (M((3) = — W g^ ) ■ 
We call M s {j3) the magnetization derived from Hamiltonian and Mg{0) the magnetization derived from QJ. From 
eqs. (|B2|) . we also obtain 

M,(J,A,h,D;P) = ys(S + 1) M s ( J, A, VS(5 + l)h, D; - j (7a) 

' M J ' ,A. " iflV (7b, 



v/5(5 + l) + ' ^(STl) , 5(S + l) ; 

From eqs. Q we obtain that Ms is a homogeneous function of degree 1 at all temperatures. This condition is satisfied 
by the /3-expansion of the magnetization, derived from eq. i|B3(l . 

Fig. Jl^) shows the magnetization M s as a function of h/J at J (3 = 0.4. We choose the same values of constants 
in the Hamiltonian J5J as in Fig. A — 1 and D/J = —0.5. Comparison of Figs. JSh) and J3Jd) shows a closer 
similarity of quantum and classical magnetization curves (even for low values of spin such as S = 3/2) than that of 
magnetic susceptibility curves, as far as high temperatures are concerned. 

The correlation function of spin z-componcnts between first neig hbors is written as (S?S? +1 ) S = for a given 

spin s. Its percental variation (when s varies by half-integer steps), shown in Fig. ©, is defined as 



: Iqzqz \ _ ^+1)0+1/2) - (S-S- +1 ) s 

°s{bi b i+1 ) = (S Z S Z +1 ) S ' ( ' 

For S ^> 1, the function S s (Sf has an expansion in S^ 1 whose leading term is S~ 3 , and its coefficient depends 
on the temperature. In Fig. JSJ) we take A = 1, h/J = 0.3 and D/J = —0.5. In the very high temperature of 
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FIG. 3: (a) The magnetic susceptibility per site Xs(JP) obtained from expansion l|B3fl for various values of spin, including the 
classical limit (solid line). We take A = 1, h/J — 0.3 and D/J — —0.5. (b) Relative percental error between Xs and Xoo, for 
several values of S. 



J/3 — 0.2, only the correlation functions for S = 1/2 and S = 1 differ from the classical curve by a difference smaller 
than 0.6%. As we lower the temperature, spins of higher values have correlation functions between first neighbors 
that differ slightly from the classical ones, but even at J/3 ~ 0.8, the correlation function for S — 2 approximates the 
classical result with an error smaller than 1%. 



III. EXTENSION OF THE THERMODYNAMICS OF THE QUANTUM SPIN-s XXZ CHAIN TO 

LOWER TEMPERATURES 

Although there are only seven terms in the high-temperature expansion of the Helmholtz free energy of the spin-s 
XXZ model (see appendix JbJ), different Pade approximants permit us to enhance our analytic results and extend 
them to lower temperatures. In this section, we shall assume a quantum spin chain with unitary norm. 

Bernu and Misguich [14| presented an approach to "interpolate" the high- and low-temperature behaviors of the 
specific heat per site, provided that the ground state energy per site of the chain is known. We refer the reader to 
their work for details on the method. 

The XXZ models (with S ^ 1/2), at the isotropic points (A ± 1), are among the most studied ones. For this 
reason, in this section we calculate the Pade approximants at those points, for the sake of comparison of our results 
to the literature. 

At T ~ the behavior of the specific heat per site is distinct for gap and gapless models. Ferromagnetic and half- 
odd integer spin chains do not have excitation energy gaps and the specific heat per site around T = is C s ~ T p / q , 
where p and q are integers|l4j. By contrast, the integer antiferromagnetic spin chains have an energy gap between 
the ground state energy and the lowest exited state and this thermodynamic function, in the region of T ~ 0, is 

C s - e--^T a ^. As the spin increases, this energy gap decreases rapidly; it only vanishes at the classical limit 
(S — > oo), though. 

In Ref. 20], Yamamoto studied numerically the thermodynamics of the S = 2 XXZ model with A = 1, D = and 
h = 0. The ground state energy per site of this model for the non- unitary spin chain was obtained in Ref. [2lJ; the 
ground state energy per site of the unitary spin chain can then be easily obtained, being equal to e /J = —0.79354. 
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FIG. 4: (a) Magnetization M s (h/ J) at J/3 = 0.4, calculated from the /3-expansion of Wa(/3), for various values of spin, including 
the classical model (S — > oo). We take A = 1 and D/J = —0.5. (b) Relative percental error of M s , for several values of S, with 
respect to Moo- 



This model exhibits an energy gap E gap /J « 0.0593 21]. At very low temperature it is expected that the specific 
heat vanish exponentially (C s ~ exp(— E gap /kT)) as T — > 0. On the other hand, the s pin wave theory gives a linear 
behavior for the specific heat over a finite interval of intermediate low temperatures |22). In order to enhance our 
/3-expansion of Ca(T) to lower temperatures, we simply use C2 ~ T. The entropy at low temperatures behaves as 
5(e) = (e — eo) 1 / 2 , where e is the ground state energy per site, and it has an essential singularity when e — ► cq. 
Following Ref. [0J, we construct an auxiliary function Q(e) = (5(e)) 2 , which is analytic in the interval [eo, 0], and use 
the Pade approximant to fit G(e), so the entropy becomes 5(e) — Q(e) 1 ^ 2 . By Q(e) we denote the Pade approximant 
realized over Q{e). The specific heat per site of the model is obtained using the relation C-z{e) = —S'(e) 2 /S"(e), 
which can be plotted in parametric form {T(e), 02(e)} (for details see Ref. [lj])- 

In Fig. (top) we compare all "acceptable" Pade approximants P n , m (i-e., all approximants so that n + m = 7 
and that do not possess spurious singularities) and their average P avr , to the numerical results obtained by the 
quantum Monte Carlo method (QMC)[2(j and the /3-expansion of Cq.(T). Fig. (bottom) displays the percental 
deviation of each Pade approximant P TOi „ with respect to P avr . We define AC2(%) = (P m ,n — Pavr)/ Pavr x 100%. 
The high temperature expansion satisfactorily describes C-ziT) up to kT/J ~ 0.9. The solid line corresponds to the 
percental deviation of P avr with respect to the QMC prediction. Although we have only seven terms in expansion 
(|B3ll of the Helmholtz free energy, the percental relative difference of the Pade approximants is smaller than 2, 5% for 
kT/J » 0.25, allowing good precision in calculating the temperature where C%(T) is maximum. 

Figs. and @ show the specific heat for several values of spin for the antiferromagnetic ( A = 1) and ferromagnetic 
(A = —1) cases, respectively, allowing us to observe how the specific heat per site approaches its classical limit as S 
varies. In order to show how the Pade approximants enhance the high temperature solution of the specific heat per 
site, Fig. also displays the /3-expansion of each C S (T). Following the spin wave theory, we assume that the specific 
heat behaves as C s ~ T for the antiferromagnetic 22] case, and as C s ~ T 1 / 2 in the ferromagnetic [2^| case. We fit 
the numerical values of the ground state energy per site calculated in Ref. |2l|, for the antiferromagnetic case with 
5 = 1,2 and 3, and extrapolate them to the integer antiferromagnetic spin S chain. Our approximate result for this 
ground state energy per site is e a /J w -(1 + 0.3641/S' + 0.029/5 2 + 0.0086/5 3 )5/(5 + 1). For the ferromagnetic case 
at the Heisenberg point (A ± 1) with h — and D — 0, the ground state energy simply becomes eo/J — —S/ (5 + 1). 

Fig. shows the specific heat as a function of temperature for A = 1, with D = and h = 0. The curves 
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FIG. 5: Percental variation of the correlation function S s {Si Sf+i) as a function of S, for several values of J (3 in the high 
temperature region. In all plots we take A = 1, h/J = 0.3 and D — —0.5. 



A = 1. D = and h = 




kT/J 



FIG. 6: (Top) The Pade approximant of the specific heat as a function of kT/J, for S = 2 with A = 1, D — and /t = 0. The 
thicker continuous line stands for the /3-expansion of C2(T). (Bottom) The AC2(T)100% represents the percental deviation of 
all acceptable (non-singular) Pade approximants and the QMC results]^] with respect to the average of Pade approximants. 



correspond to the P aV r of this thermodynamic function for the following set of spin values, as well as their respective 
/3-expansions: S — {3, 4, 5, 10, 20}. The exact classical curve corresponds to Fisher's solution^(|. From our numerical 
results, we verify that even the Pade approximants for S — 3 deviates less than 2% from the classical result (S — > oo) 
up to kT/J - 0.64, while for S = 5 this deviation is about 2% only at kT/J - 0.42. 

Except for the classical limit (S — > oo), to the best of our knowledge there are no known results of the specific heat 
per site for S > 3 for the whole range of temperature. In order to assure the validity of our results for the specific heat 
at lower temperatures, for distinct S as presented in Fig. Q, we plotted in Fig. [SI the relative percental difference 
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A = 1. h = and D = 




kT/J 



FIG. 7: The specific heat versus temperature for several values of spin at the isotropic point A = 1 with h = and D — 0. For 
each value of spin, the curve is obtained by taking the average of all "acceptable" Pade approximants (i.e., with no spurious 
singularities), whereas in the classical limit corresponds to Fisher's solution. 




FIG. 8: The relative percental difference of all "acceptable" Pade's approximants to the average Pade P aV r, for the specific heat 
per site for the antiferromagnetic case (A = 1 with h = and D = 0), with S — 4, 6, 8 and 20. Approximants with spurious 
singularities have been excluded. 



between the Pade approximants without spurious poles and the average Pade for several values of spin (5 = 4, 6, 8 
and 20). For all those values of spin, it can be seen that deviations get larger than 2% for kT/J < 0.25, whereas their 
respective /3-expansion are good approximations only up to kT/J ~ 0.8. 

In Fig. (J2J we verify that the curves of the specific heat for the antiferromagnetic case are close to the classical 
curve up to low temperatures, even for S — 3. Comparing its behavior to that of the ferromagnetic case, we plot in 
Fig. (0 the specific heat for several values of spin, taking A = — 1, h = and D = —0.1. We assume that the specific 
heat for the ferromagnetic case and small values of D still behaves as Cl ~ T 1 / 2 and the ground state energy per 
site is eo/J = —S/(S + 1) — DS/(S + 1) — D/3. The first term in the ground state energy corresponds to the case 
D = 0; the second one is the contribution of the single-ion anisotropy to the ground state energy; and the last term 
corresponds to its global shift. In order to simplify our calculations, we impose that the average energy vanishes at 
(3 = (otherwise it would be equal to D/3, that is, the value obtained from our high temperature expansion). 

As in the antiferromagnetic case, we do not know of results that would allow comparison regarding the dependence 
of the specific heat on the temperature in the interval [0, oo]. In order to carry out a first check of the curves in 
picture we proceed in a similar way as before when testing the antiferromagnetic curves. Fig. (|1C)|> shows the 
relative percental differences of the Pade's approximants and the respective /3-expansion of C S (T) to the average Pade 
for S = 1,2,3 and 5. We see that the high temperature series describe well the specific heat of the ferromagnetic 
chains for kT/J > 1. From theses plots we verify that the extension of the high-temperature results to the region of 
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FIG. 9: The specific heat per site versus temperature (kT/J) for several values of spin in the ferromagnetic case. We take: 
A = —1 with h — and D = —0.1. Each curve is obtained by plotting its P aV r, while the classical limit is obtained numerically. 




0,5 1 0,5 

kT/J kT/J 



FIG. 10: The relative percental difference of the Pade's approximants to the average Pade P avr of the specific heat per site for 
the antiferromagnetic case (A = 1 with h = and D — 0), for S = 4, 6, 8 and 20. Approximants with spurious singularities 
have been excluded. The continuous line in each picture represents the percental difference between the /3-expansion of C S {T) 

and Pavr - 



lower temperatures for the ferromagnetic case is worse than that of the antiferromagnetic case, even though integer- 
spin ferromagnetic materials are gapless and the behavior of their specific heat about T = is well described by the 
spin- wave approach. 

It is also possible to write the magnetic susceptibility per site as a function of the ground state energy per site 
e. Its low temperature behavior can be inferred from the extrapolation in a similar fashion as that of the specific 
heat per site. Unfortunately, convergence is not as much as satisfactory. By applying the method of Ref. 0] to the 
magnetic susceptibility, a large number of singularities are found in the Pade approximants to this thermodynamic 
function. For this reason we use the Dlog-Pade approximant|lq to extend the high temperature expansion of the 
antiferromagnetic susceptibility, incorporating the low-temperature information in a simple way. We follow Ref. [l5| 
to calculate the Dlog-Pade approximant of the magnetic susceptibility per site we refer the reader to this article 

for further details on the method. From Ref. [l5j, we verify that the Dlog-Pade approximants are independent of 
the leading term coefficient of the magnetic susceptibility at low temperatures, if it vanishes at T — > 0. On the other 
hand, if it is non-null at T — 0, we need to know the exact coefficient of the leading term of the magnetic susceptibility 
for T ~ in order to calculate the Pade approximants at low temperatures. This also occurs[24[ for S = 1/2. 

Using the QMC method, Yamamoto 20] obtained the magnetic susceptibility of the antiferromagnetic 5 = 2 XXZ 
model (A = 1) in the whole interval of temperature of a periodic chain with 96 sites. In Fig. (|llfl we compare the 
/3-expansion of XiiT) and its Dlog-Pade approximants for S — 2, with A = 1, D = and h = with Yamamoto's 
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FIG. 11: The magnetic susceptibility for the antiferromagnetic case A = 1, D — and h = 0. a) Dlog-Pade approximants of 
x(T) versus T/J are compared to QMC results |2C||. b) the average Dlog-Pade approximant of x{T) as a function of T/J for 
various spin values. 



numerical results. Our extension to lower temperatures assumes that the behavior of the magnetic susceptibility 
around T = be as T a exp(—E gap /kT) with a f» —1/4; the latter value yields the best fit to the QMC results p0|. 
The bottom of Fig. (|11|) shows the percental difference between the QMC and the average Dlog-Pade approximant 
to the magnetic susceptibility for S = 2; it can be seen that A\2(%) < 2% up to kT/J w 0.2. From Fig. EKt°p) we 
verify that our high-temperature expansion of X2{T) fits well the numerical result only up to kT/J ~ 0.7. 




FIG. 12: The percental difference between the magnetic susceptibility of the D-log Pade approximants and P aV r for the 
antiferromagnetic model (D — 0), with S = 1 and 5 = 3, and no external magnetic fields. 



Fig. I|12f) shows the deviation from the average P avr of the Dlog-Pade approximants to the magnetic susceptibility, 
for the antiferromagnetic Heisenberg model (A = 1 and D = 0) for for S = 1 and 5 = 3, showing the enhanced 
extention of our high-temperature results to lower temperatures. From Ref. 12 q we obtain the behavior of x around 
T = (a = 0.5); for S — 3 we assume a w 0. 

The behavior of the magnetic susceptibility at low temperature for the ferromagnetic Heisenberg model (A = — 1 
and D = 0) was obtained by Takahashi 23] using the modified spin wave theory for arbitrary spin-S. The magnetic 
susceptibility for the unitary spin chain for low temperatures, in units of k, is 



X.(T) = 2( 3 f I )V 



- - aT 1 ' 2 
3 



G(T 3 / 2 )), 



(9) 
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where a — y ^-Ji and ((x) is the Riemann zeta function. When the Dlog-Pade approximant method[lj| is applied, 

not even the coefficient of the leading term of eq.J^J) is used, but only its power. Therefore, a better approximation will 
be achieved if we take into account all the three terms in (PJ; we do so by applying the two-point Pade approximant |lfij 
method. 

We verify from expansion © that Xs(T) is singular at T = 0. Expansion l(B3() gives the /3-expansion of Xs(P) up 
to order /3 6 . The first terms of Xs(/3) for A = —1, D = and h = are Xs(/3) ~ /3/3 + 2/3 2 /9. We define the auxiliary 
function V s (0) = Xs pP — jp- This auxiliary function is regular at T — and T — > oo (/3 = 0); its only problems are 
that it has non- integer powers of /? and it does not have a Taylor expansion about T = 0. In order to circumvent this 
drawback, we proceed as in Ref. |15| and apply the transformation 



u = /3 1/2 /(l + /3 1/2 ) ^ /9 = u 2 /(l " «) 2 , (10) 

where the parameter f3 € [0, oo) is mapped onto u S [0, 1]. The expansions of V, around T = (u = 1) and /? = 
(li = 0), are polynomials in u and therefore we can write its two-point Pade's approximants, connecting the expansions 
of V s at u — 1 and u = 0. 



& = 0, -D = 0, A = -1 and S = 1 




FIG. 13: (Top) The function (kT/J) 2 xi(T) versus fcT/J of the isotropic ferromagnetic 5 = 1 model (A = -1 and D = 0) in 
the absence of external magnetic fields, for P avr - (Bottom) The relative percental difference of non-singular Pade approximants 

to Pavr- 



Fig. I|l 3(1 shows the average Pade of the magnetic susceptibility for the isotropic ferromagnetic S — 1 model, in the 
absence of an external magnetic field, obtained from the two-point Pade method and using the transformation 1(10(1 . 
For all Pade approximants, the relative difference to P avr is smaller than 0.3% in the whole interval of temperature. 

Finally, Fig. ((14(1 shows the relative percental difference of various Pade approximants to the P aV r, in the case of 
the ferromagnetic isotropic chain model (with D — 0) and in the absence of an external magnetic field. The difference 
increases with S; however, for S = 5 the relative difference is smaller than 2% up to kT/J w 0.3. 



IV. CONCLUSIONS 



In this report we present the high temperature expansion of the Helmholtz free energy of the XXZ chain, in the 
thermodynamic limit, for an arbitrary spin S up to order /3 6 . This expansion is analytic in S and in (3 and each 
coefficient of j3 n (n = — 1, 0, 1, • • • 6) is exact. The model includes anisotropy in the z-direction, single- ion anisotropy 
and an external magnetic field. In order to obtain the expansion of the Helmholtz free energy l(B3(l we apply the 
method developed in Ref. [Iflj ; from this expansion we are able to obtain the thermodynamic quantities of the model 
for arbitrary spin values. The expansion l(B3(l recovers the /3-expansion of the Helmholtz free energy presented in Ref. 
for several values of spin, and it is an extension to the results of Ref. [l^, where we include the anisotropy effects. 

The series in S allows us to obtain the classical limit of the XXZ model, but taking into account the quantum 
nature of the spins. For the particular case h = and D = we recover Joyce's result for the classical limit of this 
model[Tll|. We explicitly show in expansion Q| that only the presence of an external magnetic field distinguishes the 
classical ferromagnetic and antiferromagnetic phases in the regime of high temperatures. 
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FIG. 14: The relative percental difference between the Pade approximants to the magnetic susceptibility (ferromagnetic case) 
and P avr for S = 3/2, 2, 3 and 4. Here, A = -1, D = and h = 0. 



We interpolate the /3-expansion with its low temperature series using the Pade approximant of 1|B3|) to extrapolate 
the properties of the specific heat to finite temperatures using the method developed in Ref. |14J . Only a few terms of 
the high temperature series are needed to yield good agreement: for low spin values the deviation of all non-singular 
Pade approximants respect to P aV r fluctuates around 5%. 

For the antiferromagnetic magnetic susceptibility we use the Dlog-Pade^^j approximants to investigate the finite 
temperature properties of the spin-s quantum XXZ chain. In order to plot the magnetic susceptibility Xs(T) as a 
function of temperature for different spin values we take the average of all non-singular Dlog-Pade approximants, since 
their deviation from the corresponding P avr is less than 2%. The ferromagnetic case is studied by applying two-point 
like Pade approximants, since we have the first three terms of the low temperature series expansion 23J . The magnetic 
susceptibility is plotted for different spins, in the same way as for the previous case. 

Obviously, if we knew more low temperature information about any physical quantities, it could be also interpolated 
with our high temperature expansions, by adapting different approximants available in the literature. Our main 
limitation regarding the study of finite temperature behavior of other parameters of the Hamiltonian J2J) is the lack 
of related low temperature information. 
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APPENDIX A: NORMALIZED TRACES OF (S z ) 1 

We have that ti((S z ) n ) = if n is odd. For the purposes of this appendix, n is taken as an even integer. In order 
to calculate ti((S z ) n ) we use the result of Ref.flTj that permits us to write this trace as a function of the spin S, 

S n+l 

tr((ST)= E 3 n = Y, C ^ Sr > ( A1 ) 

j=-S r=0 

where S = 1/2, 1, 3/2, • • • , C r . n = ^ Sfe=o (r)^ fe >«+ 1 ' so ^ na ^ ^fc,«' s are the coefficients of Bernoulli's polynomial 

B n (x). They are defined bypjfl 



B n {x) 



k=Q 



(A2) 
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On the other hand, we prefer to rewrite ea. ljAip as a function of the square of the norm of the spin S (that is, 
5(5 + 1)), which is a constant of motion, 

S n/2 

£ j n = (25 + 1) A r ,„(S(S + (A3) 

j=-S r=0 

We remind that 25+1 is equal to the dimension of the Hubert space at each site, that is, tr;(l) = 25+1, i = 1, 2, • • • , TV. 
Comparing eas. (|Al(l and (|A3|) . and after some algebraic manipulation, we obtain each coefficient A r , n of ea . (|A3(I as 
a combination of the Bernoulli numbers|l7j| B n+ i-j, 

In eq. (J2J) we rewrite Hamiltonian (JJJ in terms of the rescaled spin operator s. To obtain the results l|B4|) - (|B10|) 
we need to calculate the powers of the operator s z . From the results (|A1B - (IA4|) we are able to write the normalized 
traces of s z , that is ({s z ) 2n ) just as a polynomial of y = s ^ +1 ) of degree n/2 — 1, which reads 

n/2 

((sT) = E A »/2-^y r - (A5) 

r=0 

In terms of the normalized matrix representation, the commutation relation becomes 

[s x , s y ] = —= tsZ =iJys z , (A6) 
1 J y/S(S + l) V V ' 

which is useful to evaluate the normalized traces. 

APPENDIX B: THE /3-EXPANSION OF THE FREE ENERGY FOR THE SPIN-s XXZ CHAIN 

We call W s (/3) the Helmholtz free energy of the Hamiltonian @, that drives the dynamics of the quantum spin 
chain with unitary norm. Let Ws(/3) be the Helmholtz free energy of the Hamiltonian 

JV 

H s = ]T J (Si, S i+1 ) A - hSf + D(S* ) 2 , (Bl) 

where the norm of the spin-S is \J 5(5 + 1) and 5 = 1/2, 1, 3/2, • • • . In Ref. we calculated the high temperature 
expansion of the Helmholtz free energy of the Hamiltonian (|Bl|l up to 5 = 4. The relation between the free energies 
derived from Hamiltonians and (jBlfl is 

W,(J,A,/l,D;/3) = ' + - Ws U A, </S(S + T)h, P; s{; ? + ) (B2a) 

= Ws (sisTT^Tmrrr ■ (B2b) 

In eqs. (|B2|) . the parameters of the respective Hamiltonians are shown explicitly. Equating eqs. i|B2a(> and l|B2b|) . we 
verify that Ws is an homogeneous function of degree 1. 

In this appendix we present the high temperature expansion of the Helmholtz free energy W s (/3), up to order 
/3 6 (n = 7), for arbitrary value of s(s — 5). This expansion can be written in terms of y = gr^m > that is, 



W.03) = - M2g + 1) + £ + G(/3 7 ). 

The coefficients w r (y), for r = 0, 1, • • • , 6, are shown below. For the sake of simplicity, we define x = (3y 

My) = f 

D 2 x h 2 J 2 + A J 



wi(y) 



h? _ J± _ J 2 A 2 
90 6 9 18 



,„ fVl -( j2D + D'i^V-*) _ J 2 A 2 D _ ^LDn _ J 3 A y , JAh 2 
w 2\yj \ 135 T 5670 135 90 / 36 ' 9 

... (..\ =( 2JAh 2 D _ h 2 D 2 (15 y-4) _ J 2 D 2 (27 y+4) D 4 (105 y 2 -54 y-8) _ 8 J 2 A 2 D 2 (2y-l) x 
"JZyy) \ 135 3780 18900 113400 4725 > 

h 4 (y+2) _ J 4 (-128^+12+33 y 2 ) , J 2 A 2 h 2 (3y-14) _ J 4 A 4 ( — 16 y+14+y 2 ) 
' 360 16200 270 5400 
J 2 h 2 (9 y+8) _ J 4 A 2 (-32y-72+27 y 2 ) 
+ 540 8100 



/ n _, J 2 A 2 h 2 D (16 y— 43) , 16 J A h 2 D 2 (2y-l) _ J 4 A 4 D (y 2 -12y+4) _ h 2 D 3 (105 y 2 -54 y-8) 
w i{y> —( 4725 + 4725 11340 56700 

_ J 3 AT) 2 y (9 y — 2) D 5 (15 y-4) (105 y 2 -48 y-16) _ J 4 A 2 D (3 y+4) (y-4) J 2 h 2 D (27 y+4) 

8100 + 3742200 22680 ' 18900 

, J 4 D (5 y 2 —32 y — 8) , J 2 D 3 (3 y 2 -60y-32) h 4 D (5 y+8) _ J 2 A 2 D 3 (75 y 2 -63 y+4) n 
22680 ' 170100 ~r 7560 42525 > X 

_ JAh 4 (y+2) J 3 Ah 2 (9 y+8) (y-3) , J 3 A 3 h 2 ( — 48 y+92+3 y 2 ) _ J 5 A 3 y ( — 16 y— 4+3 y 2 ) 
135 ~r 1350 ~T 4 050 6480 

J 5 Ay ( — 16 y— 4+3 y 2 ) 
3240 



/ x _ / J £> 4 (65223 y 3 +34020 y 2 -5904 y- 16064) 
UMJ/J — ( 785862000 + 8505 1 40500 

_ J 4 D 2 (4287y 3 -14910y 2 +2404y+3184) J 3 A h 2 D (39 y 2 -62 y-40) _ 2 .7 A h 4 D (16 y+27) 

23814000 ' 28350 14175 

, 2, J Ah 2 D 3 (75 y 2 -63 y+4) _ J 3 A D 3 y (213 y 2 -129 y+20) . J 2 h 2 D 2 (123 y 2 +32 y+32) 
+ 42525 425250 + 113400 

, D 6 (7182945 y 4 -6404130 y 3 +569016 y 2 +443424 y+2944) _ h 2 D 4 (15 y-4) (105 y 2 -48 y-16) 
+ 30648618000 1496880 

_ J 4 A 4 D 2 (111 y 3 -1155y 2 +888 y-16) , h 4 D 2 (105 y 2 +74 y-72) J 2 A 2 h 2 D 2 (15y 2 -36y+] 

1190700 ' 226800 + 5670 

_ J 4 A 2 D 2 (369 y 3 -546 y 2 -751 y-116) _ J 2 A 2 D 4 (204885 y 3 -224028 y 2 +39888 y +8768) 

1488375 196465500 
_ J 5 A 3 Dy 2 (4y-27) x . J 4 A 4 h 2 (-768 y 2 -2872+2748 y+45 y 3 ) _ h 6 (12y+16+3y 2 ) 

40500 > X + 340200 45360 

_ J 2 h 4 (632 y+432+333 y 2 ) _ J 6 A 2 (101120+203200 y+83481 y 2 -289584 y 3 +42174 y 4 ) 

226800 142884000 
_ J 6 (-33760- 1920 y+239418y 2 -127872 y 3 + 14067y 4 ) 



214326000 56700 
J 6 A 4 (-16960+50720 y +30846 y 2 -33624 y 3 +4239 y 4 ) , J 4 fr 2 (-108 y 2 -32-160 y+27 y 3 ) 

71442000 + 45360 

J 6 A 6 (1712-15408 y+14688 y 2 -2700 y 3 + 135 y 4 ) , J 4 A 2 fr 2 (-813 y 2 +848+383 y+135 y 3 ) 
42865200 + 85050 
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/ \ _ _ 2 .7 A h 4 D 2 x 2 , / J A 2 D (1780695 if -2280645 y 3 +677964 y 2 +89296 t^-25024) 
w 6\y) — 4725 "r I 2554051500 

,7 2 /i 2 D 3 (262035 y 3 -85176 y 2 + 12576 y- 16064) _ J 6 A 6 _D (45 iy 4 -720 »; 3 +2988 y 2 -1608 y-208) 
' 392931000 15309000 

_ J 2 h 4 D (31 y 2 +40 y+16) _ J 4 A 2 D 3 (25002 y 4 -47646 t( 3 -23907 i; 2 + 15308 j/+10464) 

75600 147349125 
_ J B A 4 D (729 y 4 -5328 g/ 3 +4488 y 2 +7784 y+1248) , ,/ 4 fr 2 D (1935 y 3 - 1974 1; 2 +2404 y+3184) 

30618000 ' 23814000 

, .7 6 A 2 D (108 y 4 -684;; 3 +201y 2 +524y+624) , J 2 A 2 h 2 D 3 (40977 y 3 -91791 iy 2 +58890 y-7240) 
' 15309000 + 19646550 

J 3 Ah 2 D 2 (7425 y 3 -11907 a 2 -2928 y-928) , J 4 A 2 ft 2 73 (972 a 3 -4851 y 2 + 1567 jy+2084) 
' 5953500 + 1786050 

_ J 2 A 2 h 4 D (35 y 2 -141 y-332) , J Ah 2 D 4 (204885 y 3 -224028 y 2 +39888 iy+8768) 

56700 ' 98232750 

_ J 5 A 3 D 2 y (1422 y 3 -8253 y 2 +3697 y-208) _ J 5 A D 2 y ( 1455 y 3 -5418 y 2 -428 y+376) 
11907000 11907000 
J 4 A 4 h 2 D (1665t( 3 -22113i; 2 +53808i/-35212) /i 4 D 3 (2625 j/ 3 +300 j; 2 - 1928 1>+224) 
' 17860500 7484400 

, J 3 A 3 /t 2 D 2 (333 a 3 -3465y 2 +4596 y-1504) J 6 D (603 jj 4 -5400 y 3 + 10062 y 2 +3520 y-416) 
' 893025 + 30618000 

_ J 4 A 4 D 3 (2y-l) (8925 y 3 -79200 j; 2 +49076 y+8912) _ J 3 A D 4 y (4413 y 3 -4214 y 2 +724 y+48) 

196465500 11907000 
_ J 2 7J 5 (150639 y 4 +58088 y 3 -28064 y 2 -67200 y- 17664) _ fr 6 D (7 y 2 +22 y+24) 
3405402000 226800 
J 4 J 3 (6201 y 4 +346320y 3 -363672 y 2 -223664^-35136) _ J A h 4 D 2 y 2 
2357586000 567 
, D 7 (15 y -4) (315315 y 4 -238140 y 3 + 1896 i; 2 + 17664 jy+2944) 
+ 30648618000 

_ h 2 D 5 (7182945 y 4 -6404130 i< 3 +569016 y 2 +443424 jy+2944) x , J Ah 6 (44 y+54+11 y 2 ) 

10216206000 ) X 28350 

_ J 7 A y (1600+6432 iy+51631 a 2 -24768 y 3 +2556 y 4 ) _ J 3 Ah 4 ,(^284 1; 2 -640-834 a+99 y 3 ) 

40824000 113400 
_ J 7 A 5 y (-400-3344 y+9228 y 2 -2964 y 3 +243 y 4 ) _ J 3 A 3 h 4 (-228 j; 2 + 1488+328 y+15 y 3 ) 
20412000 85050 
,7 5 A 3 h 2 (-33440+1080 t;+51411y 2 -25974t; 3 +2754y 4 ) J 5 A h 2 (10112+30512 t;+4629 y 2 -20844 y 3 +3159 y 4 ) 
' 7144200 + 7144200 

,7 5 A 5 h 2 (20528-31872 y+16452 y 2 -2700 y 3 + 135 y 4 ) J 7 A 3 y (3200+19808 a +14719y 2 -12912 a 3 + 1584 y 4 ) mi n"! 

' 7144200 20412000 \ OW ) 

The coefficients J, D and ft in l|B10)l are the constants in Hamiltonian 

Writing the expansion l|B3|) in terms of the variable y greatly simplifies the calculation of its classical limit (y — > 0). 
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